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Abstract 

In electrical engineering metamaterials have been developed that offer un- 
precedented control over electromagnetic fields. Here we show that general 
relativity lends the theoretical tools for designing devices made of such ver- 
satile materials. Given a desired device function, the theory describes the 
electromagnetic properties that turn this function into fact. We consider 
media that facilitate space-time transformations and include negative refrac- 
tion. Our theory unifies the concepts operating behind the scenes of perfect 
invisibility devices, perfect lenses, the optical Aharonov-Bohm effect and elec- 
tromagnetic analogs of the event horizon, and may lead to further applications. 

PACS 04.20.-q, 02.40.-k, 42.15.Eq, 77.84.-s 
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1 Introduction 



Modern metamaterials offer remarkable control over electromagnetic fields |1. with 
applications ranging from future invisibility devices PP-[Z] to existing perfect lenses 
[5]-|14|. Imagine there were no practical limits on the electromagnetic properties of 
materials. Given a desired function, how do we find the design of the device that 
turns this function into fact? In this paper, we show that general relativity provides 
clear recipes for calculating the required material properties. The practical use of 
general relativity in electrical engineering may seem surprising: relativity has been 
associated with the physics of gravitation ^5] and cosmology ^H] or, in engineering 
[T7] , has been considered a complication, not a simplification. But the design concept 
described here is rooted in a simple idea with a distinguished history: according to 
Fermat's Principle |3J CHI EH] light rays follow the shortest optical paths in media; 
they are effective geodesies, and general relativity has developed the theoretical 
tools for fields in curved geometries . Of course, here we exploit only some of the 
aspects of general relativity: we use the physics of curved space, but not the physics 
[TH] that creates space-time curvature in gravity. 1 In our case, electromagnetic 
metamaterials, not masses, create effective geometries. Our theory generalizes the 
concept behind the proposed perfect invisibility devices [TJ |2] to magneto-electric or 
moving media and it incorporates negative refraction [H]- [H|. For example, using a 
simple pictorial argument we show how to design magnifying perfect lenses. Finally, 
metamaterials may be also applied for laboratory analogs of general relativity, in 
particular for artificial black holes [20]- [SH]- In this way, we unify and generalize 
a range of physical phenomena that rely on the geometry of media and the recent 
opportunities of metamaterials. 

Metamaterials are materials with designed properties that stem from structure, 
not substance; where man-made structures determine the electromagnetic prop- 
erties, structures that are smaller than the electromagnetic wavelengths involved. 
Metamaterials have a long history: mediaeval ruby glass, for example, is a meta- 
material. Ruby glass contains nano-scale gold colloids that render the glass neither 
golden nor transparent, but ruby, depending on the size and concentration of the 
gold droplets. The color originates from a resonance of the surface plasmons [31] 
on the metallic droplets. Metamaterials per se are nothing new: what is new is 
the degree of control over the structures in the material that generate the desired 
properties. For example, in the modern version of negatively-refracting ruby glass 
[32] nano-manufactured pairs of gold-pillars on a silicon substrate generate finely- 
tuned plasmon resonances where each pair acts like a designer atom with controllable 
properties. 

Our starting point is not new either: in the early 1920's Gordon [33] noticed that 
moving isotropic media appear to electromagnetic fields as certain effective space- 
time geometries. Tamm [3U [33] generalized this geometric approach to anisotropic 
media and briefly applied this theory [33] to the propagation of light in curved 
geometries. In 1960 Plebanski [35] formulated the electromagnetic effect of curved 

1 As in other analogue models of gravity [201) we use the kinematic aspects of general relativity, 
not the dynamic ones. 
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space-time or curved coordinates in concise constitutive equations. Dielectric media 
act on electromagnetic fields as geometries and geometries act as effective media. In 
2000 it was understood [37] that the dipole forces of electromagnetic fields in media 
appear as the effects of geometries as well, electromagnetic fields act as geometries on 
media, a concept used to identify the conditions for the Abraham or the Minkowski 
momentum in media jHHj- Only very recently pp-[I] meta-material implementations 
of coordinate transformations were considered as engineering tools, ideas we take 
further here: we show that general relativity lends the most natural recipes for such 
engineering applications. 



2 Theory 

2.1 Electromagnetism in curved coordinates 

A geometry is characterized by the space-time measure (211 ds 2 = J2 a /3 9a/3dx a dx 13 , 
the metric, where we denote the space-time coordinates by x a with Greek indices 
running from to 3. Latin indices indicate the spatial coordinates and run from 
1 to 3, whereas x° = ct describes time measured in spatial units with c being 
the speed of light in vacuum. The matrix g a p, the metric tensor, may vary as a 
function of the coordinates, because space-time may be curved or because curved 
coordinates are used in flat space. The determinant g of g a p measures the 4D volume 
of an infinitesimal space-time element as \J^g times the infinitesimal coordinate 
volume [T5| I39j . We denote the matrix inverse of g a p by g a @; where, as customary 
in general relativity fSl EH]- the position of the indices indicates that g a p is co- 
variant and contra-variant under coordinate transformations. For example, 
flat space in cylindrical coordinates ct, r, (p, z is described by the metric tensor 
ga/3 = diag(l, — 1, — r 2 , — 1) with matrix inverse g a @ = diag(l, — 1, — r~ 2 , — 1) and 
determinant g = —r 2 . 

As our starting point, we use the result of general relativity [231101, derived in 
Appendix A, that the free-space Maxwell equations can be written in the macro- 
scopic form 



<9D . _ „ _ dB 

ti = 

or, in Cartesian components 



V-D = p, VxH = — +j, V-B = 0, VxE = -— 



\^ =p ' 4^=°' (2) 

E^hf^ , E^f^-f* (3) 

^— ' ox? at ^-^ oxi at 

jk jk 

where e^ k is the completely antisymmetric Levi-Civita tensor in Cartesian 

components e^ k is 1 for all cyclic permutations of e 123 , —1 for all cyclic permutations 
of e 213 and otherwise. The spatial indices indicate that in this representation 
and Hi form the components of vectors that are co-variant under purely spatial 



3 



transformations, whereas D l and B % constitute contra-variant spatial vectors. The 
charge density p and the current density j are given by \/—gj and CyJ—gj 1 of the 
four-current j a In empty but possibly curved space, the electromagnetic fields 
are connected by the constitutive equations in SI units [3o] . 



In dielectric media, the E, B vectors generate electric polarizations and magneti- 
zations that constitute the D, H fields. The constitutive equations are described 
here in implicit form, but one can also express them as D, H as functions of E 
and B |17j . In general, the electric permittivity e and the magnetic permeability 
fi are symmetric matrices - space-time appears as an anisotropic medium - but, 
in empty space, e and \x are equal, as in perfect impedance matching [4*T] . The w 
vector describes a magneto-electric coupling between the magnetic and the electric 
fields. Some materials are magneto-electric, but the simplest example is a moving 
medium [T71 EI3 E2] , because a moving medium responds to the electromagnetic field 
in locally co-moving frames and Lorentz transformations from the laboratory frame 
naturally mix electric and magnetic fields [T7| W2\ . In explicit form [36J , 



the e and fi are constructed from the spatial components of the inverse metric 
tensor and from the determinant and the time-time component of the metric tensor, 
whereas the w vector is given in terms of the time-space components of the metric 
tensor. Empty space appears as an impedance-matched anisotropic magneto-electric 
or moving medium. 

2.2 Transformation media 

Since empty space appears as a medium, what happens if a medium appears as 
empty space? This is the idea behind the recent proposals for invisibility devices 
To be more precise, suppose that the medium appears as the result of 
a coordinate transformation from some fictitious space-time, say electromagnetic 
space-time, to physical space, see, for example, Fig. [T] Electromagnetic space-time 
could be flat with light traveling along straight lines, whereas to electromagnetic 
fields physical space appears to be curved, bending light. Of course, this apparent 
curvature is an illusion, the same type of illusion as straight lines appearing curved 
in curved coordinates, because in theory it is removable by the inverse coordinate 
transformation: but, in practice, one can exploit this apparent curvature to create 
illusions 

We use primes to denote the geometry and coordinates of electromagnetic space- 
time 2 and describe physical space in generalized coordinates x l with spatial metric 
7ij and determinant 7, for keeping the theory as flexible as possible. For example, 
we may wish to use cylindrical coordinates r, ip, z with spatial metric tensor 7^ = 

2 We use the opposite notation of Ref. 1 , since, we believe, it is more logical. 



D = e £E + -xH, B 



e c 




(4) 




(5) 
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Figure 1: Transformation media implement coordinate transformations. The left figure 
shows an orthogonal grid of coordinates in electromagnetic space (a slice of cylindrical 
coordinates at constant z). The right figure shows the transformed grid in physical space 
that corresponds to the invisibility device described by Eq. (fTT|) and illustrated in Fig. [2J 



diag(l, r , 1) and 7 = r . The metric 7^ should differ from the effective geometry g a @ 
generated by the medium. The transformation rules of tensors in general relativity 
give rise to a simple recipe for the construction of media that facilitate space-time 
coordinate transformations. Contra- variant tensors are transformed as [T7J 

^£^^>' (•> 

a'P' 

The transformed inverse metric serves as the building block of the dielectric tensors 
(JHJ) of the medium. If we wish to express physical space in generalized coordinates, we 
need to consider a subtlety that appears when we write the divergences in Maxwell's 
equations (0) in spatially co-variant form 



1 cy7^ _ 1 ^ d^rB* _ 



The e u and ^ tensors are naturally contra-variant with respect to the background 
geometry of physical space, but 7 differs from —g, in general. Consequently, for 
writing the medium as an active coordinate transformation, we need to multiply D 
and B by y/j in the constitutive equations (@J) and re-scale p and j accordingly, which 
is also consistent with the form of the Levi-Civita tensor in curved coordinates 
in the curls in Maxwell's equations @. However, a second subtlety arises from the 
curls: the coordinate transformation may turn a right-handed coordinate system into 
a locally left-handed one, but the curls (j3J) implicitly assume a right-handed system, 
because e^ k changes sign under coordinate transformations from right to left-handed 
systems. Consequently, we need to invert the sign of e, \i and p, j wherever this is 
the case. The transformation to a left-handed coordinate system thus corresponds 
to negative refraction occurring in what has, for other reasons , been 
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fittingly described as left-handed media. 3 Taking all this into account, we arrive at 
the simple recipe 



e 
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y/l 900 ' 



The =F sign indicates the handedness: minus for right-handed transformations and 
plus for locally left-handed ones. Starting from the inverse metric g' al3 in electro- 
magnetic space, the g af3 matrix is calculated according to the transformation rule 
(JEJ). The matrix inverse of g al3 gives the metric tensor g a p and the inverse of the 
determinant of g a @ gives g. Equations (JSJ) specify the required electromagnetic prop- 
erties that will turn the coordinate transformation into reality. In the special case of 
right-handed and purely spatial transformations our recipe agrees with the theory 
of Refs. [HUE!) see Appendix B, where however, instead of the contra- variant e and 
\i mixed tensors occur with one of the indices lowered by the spatial metric, 



Here we consider more general transformations that may mix space and time and 
that may be multivalued. As long as the matrix (JEj) is single- valued and not explicitly 
time-dependent such a device is physically allowed and stationary. 

The coordinate transformation encodes the function of the device. Outside of it, 
the electromagnetic coordinates agree with the coordinates of physical space: the 
transformation is trivial; e and fi are unity and w vanishes. Inside the device, elec- 
tromagnetic fields are controlled as prescribed by the coordinate transformation. If 
the physical coordinates enclose a hole one obtains the blueprint of an invisibility 
device [TJ- . We show in the next section that perfect lensing [E]-[I1!, the opti- 
cal Aharonov-Bohm effect [HI 122 EH EI and optical black holes [H]-[II] can be 
understood as applications of the same idea as well. 

3 Examples 

3.1 Perfect invisibility devices 

Perfect invisibility devices [TJ [2] facilitate coordinate transformations with holes in 
physical space. In this way, electromagnetic radiation is naturally guided around 
such excluded regions. Anything placed inside is hidden. Perfect invisibility devices 
PU [2] must employ anisotropic media, because the inverse scattering problem for 
waves in isotropic media has unique solutions jJHJ US]- They also require media 
where the phase velocity of light approaches infinity at the inside of the cloaking 
layer, because coordinate transformations with holes rip apart points of zero volume 
in electromagnetic space and tear them to finite volumes in physical space, unless 

3 Negative refraction in left-handed media is not related to the case of negative phase velocities in 
gravity [431 144j , because transformations to left-handed or any other coordinates have no physical 
significance there, in contrast to transformations facilitated by media. 




(9) 
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the coordinate transformation is singular. To prove this, consider the determinant 
of the e and \x tensor (jHJ). For an invisibility device PUI2] only space is transformed, 
so w vanishes, goo is unity and — g is reduced to the inverse of the determinant of 
the spatial components . We obtain from Eqs. (JHJ) and (jHJ) in three-dimensional 
space 

w^MV-W'VJ, J = %r^4' < 10 > 

in terms of the Jacobian J jSH] ■ At a point of measure zero y/—g' vanishes and so does 
det e, the product of the eigenvalues of e and \i. Consequently, at least one of the 
eigenvalues of e and fi are zero; therefore in at least one direction of the anisotropic 
medium the phase- velocity of light diverges near the inside of the cloak. The speed 
of light in media can reach high values in narrow frequency ranges or, naturally, for 
static fields. On the other hand, invisibility devices that are only perfect within the 
limits of geometrical optics jHHHIEj are not subject to such constraints |5 a . 

Figure El illustrates a cylindrical invisibility device that stretches the z axis out in 
the radial direction to a full cylinder of radius Ri, compressed within a cylindrical 
volume of radius R%, as shown in Fig. ^ This device is an invisibility cloak of 
thickness R% ~ R\ where anything placed inside the inner radius R\ is hidden. We 
obtain from the theory 

,Ri-Ri , R2 r' /(ife-ft) 2 r 2 



within R\ <r < R2 or, equivalently, < r' < R2. Clearly, close to the lining of the 
cloak at r — ■> R\ where r' — ► the speed of light in the r and z directions diverges, 
whereas in tp direction the phase velocity tends to zero. 



3.2 Perfect lenses 

In perfect invisibility devices, electromagnetic space does not cover the entire phys- 
ical space: here we show that perfect lenses correspond to multi-valued electro- 
magnetic space. Consider in Cartesian coordinates the multi-valued transformation 
x(x') illustrated in Fig. El whereas all other coordinates are not changed. In the fold 
of the function x(x') a point x' in electromagnetic space has three faithful images in 
physical space. Obviously, electromagnetic fields at one of those points are perfectly 
imaged onto the other: the device is a perfect lens [S]. This simple pictorial argu- 
ment may contribute to settling the debate on perfect lensing jHOj-ISll i n addition 
to the experimental proof for enhanced imaging in existing perfect lenses P]-|12j. 
Our theory also shows why perfect lensing requires left-handed media with negative 
e and \i [TSUHj: inside the device, i.e. inside of the x' fold, the derivative of x(x') 
becomes negative and the coordinate system changes handedness. We obtain from 
our recipe (jHJ) the compact result 

/ dx' dx dx \ 
£ = /i = diag UPa^'d^ • (12) 

If, for example, dx'/dx is —1 inside the device and +1 outside, we obtain the stan- 
dard perfect lens [E]-[I1! based on an isotropic left-handed material with e = fx = — 1 
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Figure 2: Invisibility device. The transformation medium acts as an invisibility 
cloak, guiding light around the interior of the cloak without causing any distortion. The 
device facilitates the coordinate transformation illustrated in Fig. ^ 



inside; but this is not the most general choice. One could use an anisotropic medium 
to magnify perfect images, in contrast to the existing examples |S]-[T2], by embed- 
ding the source or the image in transformation media with |dx'/dx| 7^ 1. 

3.3 Optical Aharonov-Bohm effect 

Perfect invisibility devices and perfect lenses exploit transformation media with non- 
trivial topology, with excluded regions in physical space or folds in electromagnetic 
space. Here we consider an example where physical space is multi-valued, but the 
medium is single-valued and hence physically allowed: the optical Aharonov-Bohm 
effect [2D 1221 HH1 HZj • One can demonstrate this effect with light passing through 
a water vortex [21] or with slow light in Bose-Einstein condensates [221- Note that 
the effect is related to the Aharonov-Bohm effect with surface waves jSSJ-jSZ] and 
to the gravitational Aharonov-Bohm effect |2H] ■ The optical Aharonov-Bohm effect 
is an example of a transformation medium that mixes space and time, and yet the 
medium is stationary. Consider in cylindrical coordinates the transformation 

r' z' 
ct — ct + aip' , r = — , <p — <p' , z = — (13) 

n n 

with the constants n and a. We obtain from our theory (jHJ) that the medium 
is isotropic with refractive index n and is magneto-electric with Wi = (0, a/r, 0), 
which corresponds to a fluid vortex with velocity profile u/c = w/(n 2 — 1), as we 
see comparing Eq. (J3J) with the constitutive equations of moving media ^3 H2] in 
lowest relativistic order. Normally, a moving medium Fresnel-drags light [21], but 
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Figure 3: Perfect lens. Negatively refracting perfect lenses employ transformation 
media. The top figure shows a suitable coordinate transformation from the physical x axis 
to the electromagnetic x', the lower figure illustrates the corresponding device. The inverse 
transformation from x' to x is either triple or single- valued. The triple- valued segment 
on the physical x axis corresponds to the focal region of the lens: any source point has 
two images, one inside the lens and one on the other side. Since the device facilitates 
an exact coordinate transformation, the images are perfect with a resolution below the 
normal diffraction limit: the lens is perfect [S]. In the device, the transformation changes 
right-handed into left-handed coordinates. Consequently, the medium employed here is 
left-handed, with negative refraction |14j . 



in the case of a vortex, light rays follow straight lines, because the transformation 
f!13|) changes only time, such that straight rays in electromagnetic space are mapped 
onto straight lines in physical space. Similarly, in the original Aharonov-Bohm effect 
jSHl EDI electron rays are not bent by a magnetic vortex, but they develop a phase 
slip in the direction of incidence. In our case, when the light has passed the vortex, 
the time change in the transformation (|13|) results in a phase slip of ±nak, where k 
is the wave number, depending on whether the light propagates with or against the 
current, see Fig. |3J Physical space-time is multi- valued, with a branch cut in the 
direction of incidence, resembling the infinitely sheeted Riemann surfaces around a 
logarithmic branch point [61] . but the medium is single- valued and has the simple 
physical interpretation of a moving fluid forming a vortex. 
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Figure 4: Aharonov-Bohm effect. A fluid vortex generates the optical Aharonov-Bohm 
effect described by the coordinate transformation (|1HJI . Light, incident from the right, 
is Fresnel-dragged by the moving medium: light propagating with the flow is advanced, 
whereas light propagating against the current is retarded. The wave should develop the 
phase slip shown in the left figure. However, although the transformation is exact, 
physical space-time is multi-valued here. Instead of the simple phase slip, the light turns 
out to exhibit the characteristic interference pattern illustrated in the right figure [SHI EEJ! • 



3.4 Artificial black holes 

Moving isotropic media generate the effective space-time geometry discovered by 
Gordon in 1923 [33]. Consider an effectively one-dimensional situation where the 
medium is moving in x direction and the electromagnetic field varies along the x 
axis with field vectors pointing orthogonal to x. An impedance-matched medium of 
refractive index n and velocity u is described by the inverse metric tensor [HHUHIIE] 
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(14) 
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Both the velocity u and the index n may vary. We show that this effective geometry 
is generated from empty space by the coordinate transformation 



ct' = -(t_ + t A 



x 



t 



n ±u/c 
nu ± c 



dx . 



(15) 



as long as we restrict our attention to fields varying in x direction. For this, we 
use the transformation rule ©, but from un-primed to primed space-time, and find 
that 
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lap 



diag 



n 2 (c 2 — u 2 ) n 2 (c 2 — u 2 ) 



cf — n 2 u 2 

In electromagnetic space-time, we are free to apply the recipe 
medium properties as well as in physical space. We find e' = fi' 



9 9 ' 

n z u z 



(16) 

J) for assigning the 
1 in the directions 
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orthogonal to x and w' — 0: vacuum, which proves our point. Here wave packets are 
superpositions of waves propagating in positive or negative direction as functions of 
either t + or t_. Substituting for t± the representations ([To]) as functions of ct and 
x, we obtain the electromagnetic waves in physical space. The transformation (fTB^) 
describes the relativistic addition theorem of velocities |17| [3QJ for the medium and 
light propagating in positive or negative direction. Interesting phenomena occur 
when the velocity of the medium, u, reaches the speed of light in the medium, c/n, 
which is possible in theory and perhaps also in practice. In this case, the integral 
in t± develops a logarithmic singularity. Light propagating against the current 
freezes with exponentially increasing oscillations at a horizon [S2] , see Fig. |S| This 
horizon is completely analogous to the event horizon of the black hole if the light is 
escaping from a superluminal region and to a white hole if the light is attempting 
to enter a counter-propagating superluminal medium Horizons cut physical 

space-time into distinct regions without possible communication. They correspond 
to disconnected branches of multi- valued electromagnetic space, covering it multiple 
times, in contrast to the case of perfect lensing where the folds are connected. 
Horizons are predicted to exhibit remarkable quantum effects jB2], hi particular 
Hawking radiation [OB], that are extremely difficult to observe in astronomy, but 
may possibly be demonstrated in laboratory analogs (20J- L30 . Our method suggests 
a magneto-electric analog of the event horizon, perhaps with metamaterials, if we 
interpret the effective geometry (fTIj) as generating the constitutive equations (JU) 
expressed in terms of the D and H fields as 



Here the functions n and u of the moving medium appear as parameters of a 
magneto-electric material at rest. The material establishes horizons at u = c/n 
without singular dielectric properties, in contrast to the previous proposal [23] . gen- 
eral relativity can be put to practical use in electrical engineering, but electrical 
engineering may be also used for demonstrating some elusive effects of general rela- 
tivity. 

4 Conclusions 

Perfect invisibility devices (HI2j; perfect lenses [E]-[lI], the optical Aharonov-Bohm 
effect |2U 1221 GE1 El and artificial event horizons [20]- [SHI are all examples of one 
unifying concept: electromagnetic media that facilitate coordinate transformations. 
Adopting ideas from general relativity [121-110] we developed a concise formalism 
for finding the properties of meta-materials that turn a desired function into fact. 
We extended the previously reported method [H |2] to media that act as space-time 
transformations and that may exhibit negative refraction and signs of multi-valued 
physical space. The most interesting properties of such transformation media seem 
to stem from non-trivial topologies: spaces with holes in the case of invisibility 
devices, coordinate folds for negative refraction, multi-valued physical space for the 




(17) 
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Figure 5: Event horizons. Space-time diagram showing light propagation near two 
artificial event horizons, a white-hole horizon at the left and a black-hole horizon at the 
right. The horizons could be made by a medium moving from the right side to the 
left. The places x where the medium velocity reaches the local speed of light in the 
medium are the horizons. At the black- hole horizon the medium turns from subluminal 
to superluminal velocity; at the white-hole horizon it returns to subluminal speed. The 
exponential confinement near the horizons translates into exponential red shifts of optical 
oscillations at the black hole and exponential blue shifts at the white hole. Magneto- 
electric media could mimic the moving medium generating the horizons. 



Aharonov-Bohm effect and multiple sheets of electromagnetic space-time in the case 
of artificial event horizons. Our theory can be extended in at least two directions. 
So far we assumed that electromagnetic space-time is empty, but one could easily 
incorporate media with non-trivial e' and y! here. We also assumed that the mapping 
between electromagnetic fields in real and fictitious space is exact. In practice, the 
accuracy of geometrical optics [IE] is often completely sufficient; one could further 
generalize our method to include transformations that are only exact within the 
validity range of geometrical optics. Optical conformal mapping jSJ E] is such an 
example. These more general transformation media act as local transformations of 
the dispersion relation of light waves. They all share the same spirit: electromagnetic 
media act as effective geometries. 
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Appendix A 



In this appendix we show that the free-space Maxwell equations in generally co- 
variant form are equivalent to Maxwell's equations in a material medium with con- 
stitutive equations (jlj),©. The generally covariant Maxwell equations are 



'\ftv;X] 



<9 [A F H = 0, eo^r 



-g 



(Al) 



where F^ u denotes the electromagnetic field tensor, j M is the four-current, the square 
brackets denote antisymmetrization and the semi-colon indicates covariant differen- 
tiation. We employ Einstein's summation convention over repeated indices. The 
covariant tensor F^ v contains the E and B fields in the usual manner: 
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E x 


Ey 


E z \ 


—E x 





-cB z 


CBy 


~Ey 


cB z 





-cB x 


\-E z 


-cB y 


cB x 


o ) 



(A2) 



We define a quantity H^ v with contravariant indices by 



H 



-gF^ = eaV^g^gT'Fx,, 
1 



--g^g vp H xp 



and regard as being constructed from D and H fields as follows: 



(A3) 
(A4) 
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~D y 


-D z \ 
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-H,/c 
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H z /c 
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-Hy/C 
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Then, introducing a new four-current J M = 
be written 

d [X F H = 0, d v H^ = J» 



gj^, Maxwell's equations 



(A5) 

ll can 

(A6) 



which are the electromagnetic equations in a material medium, described by the 
constitutive equations (|A2J) - (|A5|) . with free charge and current densities J M . 

To obtain relations between the vector fields D, H and E, B consider first the 
components F 0i ; from Eqs. (jA2|) . (jA4j) and (jA5|) one obtains 



E; 



{9io9jo — 9ij9oo) Dj — 



£oV-9 



■■[jkl}g 0j gikHi, 



(A7) 



where [jkl] denotes the completely antisymmetric permutation symbol ^Hj with 
[xyz] = +1. We simplify this result as follows. The identity 



g,xg Xu = s; 
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gives 



AO 



9i\9 



goxg Xi = o 



9i0 = -goo MO- 



JO 



,'0 



9oo 



9 l3 9jo ■ 



9 3 x9 M = 9 3 o9 0i + 9 3k 9 ki = 5). 



Use of Eqs. (jA8|) or (jA9|) in Eq. ()A10|) produces the two relations 



.00 y y 



ik 



9kj = 0j, g g kj 



goo 



gkogjo = S 



(A8) 

(A9) 
(A10) 

(All) 



which reveal inverse-related 3x3 matrices. In view of Eqs. ()A10|) and (jAlllh 
transvection of ()A7|) by g b yields 



A = --^-^-g^E 3 + —[ijk]gj H k 

goo goo 



(A12) 



the first of the constitutive equations (jH) with (JHJ). 

To obtain the second constitutive relation, we employ the tensors dual to 
and H^ u . This requires use of the 4D Levi-Civita tensor [15J, given by 



-p.uXp 
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[pLuXp] , [0123] 



The dual tensors *F flu and *H tiU are defined by [To] 

1 



c pvXp p 



*H, 



1 



fJ,V 



-puXp 



Xp 



pv 



-pvXp 



*pXp 



H 



p,v 



-pvXp 



* H Xp 



(A13) 

(A14) 
(A15) 



so they have components 



*ppu 



*H, 



pv 



-g 



1 


-cB x 


-CBy 


-cB z \ 


cB x 





E z 


~Ey 


CBy 


-E z 





E x 


\cB z 


Ey 


—E x 


o / 


( 


H x /c 


Hy/C 


H z /c\ 


-H x /c 





D z 


~Dy 


-H y /c 


-D z 





D x 


\-H z /c 


Dy 


-D x 


o ) 



(A16) 



(A17) 



Re-expressed in terms of the dual tensors, the constitutive equations (|A3j) - ([A4)) read 

= £ 0V ^ggpxg,p *F Ap , (A18) 
l 



*ppu 



£ox/-g 



--g»Y p *Hxp. 



(A19) 
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Writing out *H 0i using Eqs. (JA16|) - (JA18j) one finds 

2 

Hi = -"^7= (9oo9ij - 9io9jo) Bj + —f=\jkl)g j0 g ik Ei . (A20) 
Comparison of this with Eqs. (|A7|) and (|A12|) shows that 



Bi = --^-g* 3 ^ [ijk]g j0 E h , (A21) 

which is the second of the constitutive equations (j3J) with 

Clearly, several other relations between D, H, E and B are contained in (|A3jl - 
flA4|) and ()A18|) - ()A19|) . For example, to express D and H in terms of E and B, as is 
done in Eq. (fT7j). we need only take the time-space components of (|A3[) . obtaining 

A = ^ov^ foV° - <W°) ^ - eoc^L^l^V^ , (A22) 
and the required formulae are Eqs. (|A20|) and (|A22[) . 



Appendix B 

In this appendix we show that the expressions for e 1 ^ and /i*- 7 derived in Ref. [13] and 
utilized in Ref. [T] are a special case of the formalism used in this paper. Consider 
the effective e l i and // J ' resulting from a transformation of the spatial coordinates of 
a Minkowski system. From Eq. JSJ) we obtain 

e* = -V=^V, A** = -V=^V, (Bl) 

where we have allowed for general isotropic permittivity and permeability e' ^ 1, 
jj! ^ 1 in electromagnetic space. 

Rescale the spatial coordinate basis vectors <9j so that they form a (in general 
non-coordinate) basis Uj of vectors of unit length: 

Uj = —= di , g(uj, Uj) = 1 . (No summation.) (B2) 
V~9a 

Let g~ij be the components of the metric tensor in the basis u^: 

g~ij = g(uj, Uj) = Qjj . (No summation.) (B3) 

\J 9ii9jj 

We need to compute the triple product ui ■ (112 x U3) in the coordinate basis for 
this we require the 3D Levi-Civita tensor in this basis: 

e ij k = V II 9[ijk], e ijk = --^[ijk], (B4) 
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where g is both the determinant of the space-time metric and the negative deter- 
minant of the spatial metric since g 00 = 1> 9io — 0. Using Eqs. ()B2|) - ()B4|) we 
find 



Ui • (u 2 x u 3 ) 



-gijU x e 3 u 2k U3i = -giju\s J g km u 2 g in u 3 



1 iW 1 



V~9i 



1 



-913 



V-9ll922g33 \T r 9 



V~922 V~933 

[jkl}gijg2kg3i 



y/ -911922933 v 7 -^' 



V / -5 , Hfi'22fi'33 

Using Eqs. JEHJ) and (jB5j) we write (JBTJ) as 



(B5) 



1 



e % ' = -e'ui • (u 2 x u 3 )V-PiiP22^33 — : r ' 



\J 9ii9jj 



H 13 = -//'ui ■ (u 2 X U 3 )V-5'll5 , 225 , 33 



y9ii9jj 



(No summation.) (B6) 
(No summation.) (B7) 



which are the expressions derived in Ref. jlo]. Note that this form of the theory 
implicitly contains the case of negative refraction, because for a transition to a left- 
handed system the cross products change sign, but this fact was never mentioned 
nor used so far. 
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